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The low-mode averaging technique is a powerful tool for reducing large fluctuations in correlation 
functions due to low-mode eigenvalues of the Dirac operator. In this work we propose a general- 
ization to baryons and test our method on two-point correlation functions of left-handed nucleons, 
computed with quenched Neuberger fermions on a lattice with extension L = 1.5 fm. We show 
that the statistical fluctuations can be reduced and the baryon signal significantly improved. 
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1. Introduction 

It is a known fact that arbitrarily small eigenvalues of the massless Dirac operator D can occur 
[jl|, ||, ||, |], |5| and that local "bumps" of the corresponding wave functions can be reflected in 
large fluctuations of the physical observables & In particular, in the quark mass region m ~ 1 /LV, 
where V is the volume and £ is the bare quark condensate, we expect that few eigenvectors could 
give a substantial contribution to the observables; in that case the fluctuations can be reduced in 
an efficient way by applying an exact low-mode averaging procedure. By adopting the techniques 
developed in [{7|, ||[|, low-mode averaging has been applied to quenched meson correlation functions 
with Neuberger fermions, both in the e- and /^-regimes [|| §]. 

This technique can be generalized to a wider class of correlators; in particular this work is devoted 
to its application to baryon two-point functions. After clarifying the theoretical framework, we 
present a numerical study with quenched Neuberger fermions in the /^-regime. 



2. Low-mode averaging for baryonic two-point functions 



In the following we consider a lattice of volume V = L 3 T with lattice spacing a and peri- 
odic boundary conditions in all directions. We assume that fermions are discretized by using the 
Neuberger-Dirac operator D [10], which satisfies the Ginsparg-Wilson relation [11]; this ensures 
that chiral symmetry is preserved at finite lattice spacing Jl2|]. The conventions used in this work 
are the same as in [[§], to which the reader can refer for undefined notations. We adopt the neutron 
interpolating field 

N L (x) = [uf {x)Cy 5 d h L {x)\ d c L (x)e abc , (2.1) 
where q= (l — \aD\ q and C is the charge conjugation matrix. We consider the two-point function 

C N (t) = £(Tr (N l (x)N l (0)y )), (2.2) 

X 

where the trace is meant over the nucleon spinor indices. Following ref. tfjty, for each gauge 
configuration we can extract the first n low modes of the Dirac operator and express the left-left 
propagator as the sum of the light and heavy parts: 



where P+ 



S L (x,y)=P-S(x,y)P+=P- 
;i±7 5 )/2and 



1 £ ( p i (x)( P J(y)+S h (x,y) 



i=0 



+ ! 



(pi{x) 



1 



[P c \l/i + P- c DP c Yi](x). 



(2.3) 



(2.4) 



Here, P c (c = ±) is the projector on the chiral sector with no zero modes, and A,-, ifi are respec- 
tively the approximate eigenvalues and eigenvectors of A = P c D] n D m P c , where D m is the massive 
Neuberger-Dirac operator. The use of the left-handed propagator avoids the complications due to 
the contribution of the zero modes to the correlation function. The baryon correlation function in 
eq. (2.2) can be split into four contributions: 

C N (t) = C't(t) + Cf l (t)+C h ^{t) +Cf{t), (2.5) 
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where h (I) means that the heavy (light) part of the propagator is involved. The ensemble average of 
each contribution is translational invariant, and therefore its Monte Carlo variance can be reduced 
by averaging over equivalent space-time points. There are two possible contractions that contribute 
to the correlation function of eq. {12), and in the case of degenerate quark masses, the /// part can 
be written as 



<#(0 = ^E&*-* I Tr \p+ U{x)<pf \y)) T CysP-tfixWp^Cy, 

v x,y i,j=0 L 



(2.6) 



^Tr\p^ k (x)cp e k \y)y 



k=0 



£ abc £ feg 



where the translational invariance can be fully exploited without any additional computational cost. 
For the Ml contribution, the direct application of low-mode averaging would require the inversion 
of the heavy propagator on 0{n 2 ) source vectors. The strategy is then to select in a systematic 
manner only a restricted number of contributions on which space-time averaging will be applied. 
The first step in this direction is to introduce diquark vectors 



(2.7) 



where i = i\+n ii, and a,j8 are the explicit spinor indices. In the space spanned by the V we 
define the scalar product 



x,a,a,p,a' ,j3' 



(2.8) 



and we define an orthonormal basis V by diagonalizing the metric matrix {V ,L y ). The basic 
building block 



J)(x,y)=£ 1 L*(x)L*t(y) 
k=0 



can be rewritten as 



n l -\ 



S(x,y)= £ \c k \ 2 V k (x)V k Hy), 

k=0 



(2.9) 



(2.10) 



where |co| > \c\ \ > ... > \c n 2_\ \ are appropriate coefficients. The Ml contribution can be expressed 
as 



Cfit) 



1 



I^-2lk| 2 Tr 

k=0 



P.S\x,y)Cy 5 P + V kJf (y) ) Cy 5 P_V k (x)P.y 



+ 



(2.11) 



n l -\ 



I k,| 2 Tr 



k=0 



PS h (x,y)P + Tr(Cy 5 V ki! (y)P+ 



7o 



Tr 



V k (x)P^Cy 5 



and low-mode averaging is applied only on the (supposedly) dominant contributions (k = 0, ...,nj, — 
1), while the remaining terms (k = nt,, ...,n 2 — 1) are computed locally. The other contributions, Mil 
and hhh, are computed locally as well, under the assumption that their fluctuations are much milder 
than in the previous two cases. 
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Figure 1: The coefficient of eq. (2.10) as a function of k, for a given configuration at am = 0.04. 
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Figure 2: Monte Carlo history of the correlation function C„(t) at t = yo = 9a, for am = 0.04. 

3. Numerical results and discussion 

We tested our method on quenched Neuberger fermions defined on a lattice with L/a = 12, 
T/a = 32, jS = 5.8485, corresponding to a ~ 0.12 fm, L ~ 1.5 fm. We used three different quark 
masses am = 0.04,0.055,0.07 and collected N con f = 109 measurements of the baryon correlation 
function. We extracted n = 20 low modes of the Dirac operator, and chose n\, = 20 for the low-mode 
averaging on the Ml contribution. Taking into account the symmetries of the source vectors V k , this 
corresponds to computing 2«^ inversions of the propagator. In fig. [j] we show the coefficients of 
eq. ( 2.10| ) for k = 0, ...,«£, — 1 for a typical configuration at am = 0.04; they drop off considerably 
fast, which justifies the restriction k = 0, ...,ni, — I for the application of low-mode averaging. In 
any case we stress that we performed no systematic study to optimize the choice of the parameter 
rib- In fig. [2] we show the Monte Carlo history of C^{t) at a given time slice t = yo = 9a, for 
the lightest mass am = 0.04. The first figure from the left corresponds to the local computation, 
and the presence of large fluctuations can be observed. The figure in the centre is obtained by 



performing low-mode averaging on the /// part of the correlation function, as in eq. ( |2.6j ): the 
fluctuations are reduced, but nevertheless we still observe residual spikes. Finally, the third curve 
on the right corresponds to the computation where low-mode averaging is applied to both /// and hi I 
contributions: in this case the large fluctuations are sensibly suppressed. This confirms a posteriori 
that the fluctuations given by the /// and Ml contributions are the dominant ones. 
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am 


am e ff (no lma) 


am e ff (lma on ///) 


am e ff (lma on lll+hll) 


0.04 


0.69(62) 


0.76(14) 


0.78(7) 


0.055 


0.89(20) 


0.81(11) 


0.85(6) 


0.07 


0.94(11) 


0.87(8) 


0.90(5) 



Table 1: Numerical results of low-mode averaging (lma) for the effective mass in units of lattice spacing at 
t/a = 9, for different values of the quark mass am. 
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Figure 3: Effective mass for am = 0.04. 



At large time separations, we expect 



C N (t) 



A N (e- mNt -e- m ^ T -^ +A N * (e^* 1 - e-' 11 ^ 7 ^ 



(3.1) 



where m^ is the mass of the ground-state positive-parity nucleon, and m^* is the mass of the 
negative-parity state N* K In table |3| we report the numerical results for the effective mass at 
t/a = 9, for our three values of the quark mass. For the local correlator, the fluctuations are such 
that a jackknife statistical analysis of the errors cannot be trusted, and the uncertainties reported 
in the first column for am = 0.04 and 0.055 are just for completeness. For the lightest-quark 
mass, low-mode averaging applied on /// and Ml with n/, = 20 guarantees a reliable determination 
of statistical errors. A brute comparison of the errors indicates a reduction by a factor 10, to 
be compared with an increase of the computational cost by a factor 8. Going to heavier quark 
masses the low-mode averaging becomes less efficient, as expected; the same effect is foreseen by 
increasing the volume at fixed quark mass. In fig. |3| we report the results obtained for the nucleon 
effective mass, for am = 0.04. On the first figure on the left, no low-mode averaging has been 
applied, and the signal for the nucleon is very poor. In the centre we show the effect of applying 
low-mode averaging on the /// part only: for t/a >7we already observe a significant improvement 
on the signal. Finally, on the right we show the effective mass that we obtain after applying low- 



1 Notice that, because of our choice of the interpolating field, we cannot distinguish between different parities, hence 
we have both contributions of A' and N* in both time sectors. Experimentally A m = mjy* ~~ m N ~ 600 MeV, and we expect 



that the contribution of N in eq. (3.1) dominates at sufficiently large / 
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mode averaging on /// and hll contributions; here the statistical errors in the range t/a ~ 9-12 are 
under control and we can attempt to extract the nucleon mass my. It is important to stress that our 
uncertainty on the nucleon mass after low-mode averaging is still quite large with respect to the 
errors obtained without low-mode averaging but with different choices of the nucleon interpolating 



field (see e.g. [ |13| ] for a recent computation). For this reason, as a next step we intend to apply 
this method to baryon two-point functions with other interpolating fields. Moreover, the technique 
can be easily generalized to baryon three-point functions and could therefore be very useful for the 



computation of the electric neutron dipole moment [14, 15]. 
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